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Abstract 


In  this  report  we  show  how  to  construct  analytic  solutions  of  the  Klein-Gordon 
equation  in  a  semi-infinite  channel.  The  Klein-Gordon  equation  can  be  derived  from 
the  shallow  water  equations.  The  analytic  solutions  are  given  for  various  choices  of 
initial  and  boundary  conditions. 
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1  Statement  of  the  Problem 


Consider  the  shallow  water  equations  (SWEs)  in  a  semi-infinite  channel.  For  simplicity 
we  assume  that  the  channel  has  a  flat  bottom  and  that  there  is  no  advection,  although 
these  assumptions  may  be  removed  in  future  studies.  We  do  take  into  account  rotation 
(Coriolis)  effects.  A  Cartesian  coordinate  system  (x,  y )  is  introduced  such  that  the 
channel  is  parallel  to  the  x  direction,  as  shown  in  the  figure.  The  width  of  the  channel 
is  denoted  b. 
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Figure  1:  Setup  for  the  wave-guide  problem  in  a  semi-infinite  wave  guide 

The  SWEs  are  (see  [1]): 

d tu  +  gudxu  +  gvdyu  —  fv  =  —g  dxg  , 

dtv  +  gudxv  +  gvdyv  +  fu  =  —g  dyg  ,  (1) 

dtr 7  +  gudxg  +  gvdyg  +  ( h0  +  gg)  (5 xu  +  dyv)  =  0  . 

Here  t  is  time,  u(x,  y,  t )  and  v(x,  y,  t)  are  the  unknown  velocities  in  the  x  and  y  direc¬ 
tions,  ho  is  the  given  water  layer  thickness  (in  the  direction  normal  to  the  xy  plane), 
g{x,y,t)  is  the  unknown  water  elevation  above  ho,  f  is  the  Coriolis  parameter,  and  g 
is  the  gravity  acceleration.  We  use  the  following  shorthand  for  partial  derivatives 


The  parameter  g  is  1  for  the  nonlinear  SWEs,  and  is  0  for  the  linearized  SWEs  with 
vanishing  mean  flow.  We  shall  only  consider  the  latter  in  the  sequel. 

It  can  be  shown  (see  [2])  that  a  single  boundary  condition  must  be  imposed  along 
the  entire  boundary  to  obtain  a  well-posed  problem.  On  the  south  and  north  channel 
walls  T 5  and  Tjv  we  have  v  =  0  (no  normal  flow).  On  the  west  boundary  Tf  we 
prescribe  g  using  the  Dirichlet  condition  g(0,  y,t)  =  riwilht),  where  gw(y,t)  is  a  given 
function  (incoming  wave).  At  x  —>  oo  the  solution  is  known  to  be  bounded  and  not  to 
include  any  incoming  waves.  To  complete  the  statement  of  the  problem,  initial  values 
for  u,  v  and  g  are  given  at  time  t  =  0  in  the  entire  domain. 

It  is  easy  to  see  (e.g.,  [3])  that  the  system  (1)  is  equivalent  to 

dtV  ~  Cffl2g  +  f2g  =  0 


(2) 


where  Co  =  y/gfiQ.  The  boundary  conditions  are 


V(0,y,t)  =  F(y,t), 

dnv{x ,  0,  t)  =  0, 

dnr]{x,  b,t)  =  0,  ‘  ' 

lim  ri(x,  y,  t )  is  bounded 

X—>00 

and  the  initial  conditions  are 

V(x,y,0)  =  G(x,y),  () 

dty{x,  y,  0)  =  0 

The  general  solution  to  the  problem  can  be  found  by  taking  the  Fourier  sine  transform 
in  x  and  then  solve  the  resulting  PDE  in  t  and  y.  This  is  very  messy  and  will  not  be 
pursued  here.  In  the  following  sections  we  will  discuss  the  possible  choices  for  F(y,  t ) 
and  G(x,y).  One  can  also  take  a  non  vanishing  second  initial  condition. 


2  Construction  of  Solutions 


Since  the  problem  (2)-(4)  has  a  non-homogeneous  boundary  condition,  we  will  decom¬ 
pose  7] 

r)  =  v  +  w  (5) 

where  v  will  satisfy  the  PDE  with  a  homogeneous  boundary  condition  on  the  west 
boundary  and  w  satisfies  the  non-homogeneous  boundary  condition  with  a  simplified 
PDE.  We  will  choose  w  to  satisfy  the  following  problem 

8 -  CfidyW  +  f2w  =  0 
w(0,y,t)  =  F(y,t), 
dniu(x,  0,  t)  =  0, 

dnw(x,  b,  t)  =  0,  (6) 

lim  w(x,  y,  t )  is  bounded 

x — >oo 

w(x,y,  0)  =  H(y), 
dtw(x,y,0)  =  0 


The  simplification  is  the  fact  that  w  is  independent  of  x. 

Remark:  If  the  boundary  condition  at  x  =  0  is  a  function  of  y  only,  then  the  above 
equation  (6)  will  not  work.  We  will  consider  that  case  in  the  next  section. 

Problem  (6)  can  be  solved  using  separation  of  variables,  i.e.  by  assuming  that 
w  =  Y  (y)T(t). 

It  is  easy  to  see  that  the  solution  is 


where 


w  = 


OO 

E 

771=0 


Am  cos( vo(m)  t)  cos 


vo  (m) 


/2+(^)2 


(7) 

(8) 


Since  w  is  independent  of  x,  this  is  also  F(y,t).  Using  the  initial  condition 


H{y)  =  Am  cos  ( -7 -yj 

m= 0  \  0  / 

we  can  find  the  Fourier  coefficients  Am .  Let’s  choose  (for  simplicity) 

H{y )  =  cos  (j-yj  . 

This  choice  will  simplify  the  computation  of  the  Fourier  coefficients  Am, 

A  =  I  1  m  =  1 

m  [ 0  m/1 

Therefore 

F(y,t)  =  cos(^/z/0(l)  t)  cos  ^yj 

Clearly  if  we  choose  a  different  H(y),  we  get  a  different  function  F(y,t). 
w  is  then  given  by 

w  =  cos(^)  cos 

Now  we  take  v  given  by  (5)  and  substitute  in  (2),  we  have 

dfv  +  d?w  -  CqV2v  -  CqV2w  +  f2v  +  f2w  =  0 
Move  all  the  terms  with  w  to  the  right  and  use  (13),  we  have 

d?v  -  CqV2v  +  f2v  =  0 

which  is  identical  to  (2).  Now  the  boundary  conditions  become 

v(0,y,t)  =rf(0,y,t)-w(0,y,t\  =  0 

F(l l,t)  F(y,t ) 

dnv(x,  0,  t)  =  0 
dnv(b,y,t )  =  0 
lim  v(x,  y,  t )  is  bounded 

x—>oo 

The  initial  conditions  are 

v(x,  y,  0)  =  G(x,  y)  -  w(x,  y,  0) 

c°s(f  y) 

dtv(x,y,  0)  =  0 

To  solve  (14)-(16),  we  use  separation  of  variables 


(9) 

(10) 

i.e. 

(11) 

(12) 

The  solution 

(13) 

(14) 

(15) 

(16) 


v(x,y,t)  =  T(t)<t>(x,y) 


(17) 


Substituting  in  (14),  we  have  two  differential  equations 


T  +  z/T  =  0 


(18) 


with 

T(  0)  =  0 

and 

•?  v  ~  /2 
v  4>  h — ^2 — 0  =  o 

with  the  boundary  conditions 

0(0,  y)  =0 

dy<j){x ,  0)  =  0 

dy(j)(x,  b)  =  0 
lim  (j)(x,  y)  is  bounded 

>oo 


(19) 

(20) 


(21) 


To  solve  (20),  we  will  separate  the  variables  again,  assuming  (fr(x,y) 
get 


Y"  +  fiY  =  0 


X"  + 


^-/2 

C'o2 


X  =  0 


y'(o)  =  o  x(o)  =  o 

T'(6)  =  0  lim  X(x)  is  bounded 

x — >00 


The  solution  for  the  Y  equation  is 


X(x)Y(y)  to 


(22) 


(  mir\2 

*im=VTJ 

Ym(y)  =  COS 

m  =  0, 1, 2, . . . 


(23) 


In  order  for  the  X  equation  to  have  a  non  trivial  solution,  we  must  have  u  >  vq (m) 
where  v^irn)  is  given  by  (8).  In  this  case  the  solution  will  be 


s  ■  (  vV  -  Mrn) 

X  (x)  =  sm  - — - x 


C0 


Therefore  the  solution  for  (20)  is 

OO 

4>(x,y)  = 


m= 0 


'  vo  (m) 


A  f  \  •  (  \JV  ~  Mm)  ,  j 

Am[y)  sm  - — - x  \  dv 


Co 


cos 


(  mir 
\~l 


The  solution  of  the  T  equation  is 


T{t )  =  cos  {yfit) 


and  therefore 

OO 

v(x,y,t )  =  ^ 

m= 0  L' 


/  v  o  (m) 


Am(y)  cos  {s/ut)  sin  f - ~C~~ — ~  x  J  (^u 


(24) 


(25) 


(26) 


/ U7.7T  ,  ,  , 

cos  f  ■ —  y  )  (27) 


The  only  condition  left  to  satisfy  is  v(x,  y,  0)  =  G(x,  y )  —  cos  (jy)  ■  Before  we  do  that, 
let  us  transform  our  general  solution,  by  taking 


A  = 


\Jv  -  t'o(ra) 

Co 


(28) 


We  have 


(x,  y ,  t)  =  ^2  J  Bm(A)  cos  CqA2  +  uq (to)  t  \  sin  (A  x)  2C$AdA 


m= 0 


(  mn 

cos  I  —  y 


At  t  =  0,  we  have 


G{x,y) -cos(-y)  =  £ 

771= 0 


Bm(A)2(7oAsin  (Ax)  dA 


/  to7t 

cos  I  —  y 


Let  us  choose  for  simplicity 

G(x,  y)  =  g(x)  cos  (j^yj  +  cos 


then 


m= 0  CO 


-Bm(A)2C*oAsin  (Ax)  dA 


/  to7t 

cos  I  —  y 


(29) 

(30) 

(31) 

(32) 

(33) 


ff(x)  COS  l^  -yj  =  Y 
and  therefore  to  =  1  and 

r  oo 

g(x)  =  2Cq  /  ABi(A)  sin  (Ax)  dA 
Jo 

This  means  that  2CqAB\(A)  is  the  Fourier  sine  transform  of  g(x).  Let  us  choose  g(x) 
as  (the  choice  should  be  such  that  the  Fourier  sine  transform  of  this,  yield  a  convergent 
integral  in  (29)) 

g{x)  =  x 


then 


x2  +  a2 
Aa 


5l(A)  2CqA 

Now  we  substitute  this  into  (29)  to  get 

v(x,y,t)=  e~Aa  cos  ^  C%A2  +  z'o(l)  tj  sin  (Ax)  dAj  cos  y 
Combining  this  with  the  solution  for  w(x,  y,  t )  given  in  (13)  we  have 


(34) 

(35) 

(36) 


y{x,y,t)  = 


cos 


\J v o(l)  tj  +  J  e  Aa  cos  ^Y^C^A2_T^o(l)"t^  sin  (Ax)  dA  cos  y 

(37) 


G(x’y)  =  (?f^  +  0 cos  (t) 

F(y,t)  =  cos  )tj  cos  (j;yj 


(38) 


This  solution  assumes 


3  The  case  that  F(y,  t)  is  independent  of  t 


In  this  case,  one  cannot  use  (6)  because  the  solution  (13)  depends  on  time.  Instead  of 
(6),  we  should  take 

— Cg  (d2w  +  dyW^j  +  f2w  =  0 

subject  to 

w{0,y)  =  F(y), 
dnw(x,0 )  =  0, 
dnw(x,  b )  =  0, 
lim  w(x,y )  is  bounded 


(39) 


(40) 


The  solution  is  given  by 


OO  /  \ 

w(x,y )  =  y4me_("'o(m)/Go);Ecos 

m=  0  \  0  J 


(41) 


where  vq (m)  is  given  by  (8)  and 

OO  / 

F{y)  =  2^  Arn  COS  y  —  y 

If  we  choose 


m= 0 


then 


F(y)  =  cos  (j-yj 
w(x,y )  =  cos  {^V^j 


(42) 


(43) 


(44) 


This  w(x,y)  will  give  the  same  PDE  for  v.  The  only  condition  affected  is  (16)  which 
is  now 


v(x,y,0)  =  G{x )  —  e  ("o(i)/Co)a:  cos 
For  simplicity,  we  assume  that 

G(x)  =  g{x)  cos  (^yj 

and  therefore 

v(x,y,  0)  =  (jj(x)  —  cos  (^y 

This  condition  is  only  used  when  we  reach  (30)  where  we  now  have 

(  \  oo  ■  roo 

=  E  [J  Sm(A)2C'g  A  sin  (A  x) 


(45) 

(46) 

(47) 


771=0 


dA 


(  m.ir 

cos  I  —y 


Therefrore  m  =  1  and 


roo 

g(x)  -  e-toW/Cb)*  -  i  =  /  Bl{A)2C$Asm{Ax)dA 

Jo 


(48) 

(49) 


If  we  now  take 


x 


9{x)  = 


_|_  e-(^o(l)/Co)a:  _|_  i 


x2  +  a2 

then  B i(A)  is  given  by  (35)  and  v(x,y,t )  is  given  by  (36)  as  before, 
chose 

G(x,y)  =  (^2+e-^(1)/^  +  l)cosgy) 


^(y) 


(50) 

In  this  case  we 


(51) 


Acknowledgments 

The  authors  would  like  to  express  their  gratitude  for  the  support  extended  to  them 
by  the  Naval  Postgraduate  School. 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


References 


[1]  J.  Pedlosky,  Geophysical  Fluid  Dynamics ,  Springer,  New  York,  1987. 

[2]  B.  Gustafsson,  H.-O.  Kreiss  and  J.  Oliger,  Time  Dependent  Problems  and  Differ¬ 
ence  Methods,  Wiley,  New  York,  1995. 

[3]  V.  J.  van  Joolen,  Application  of  Higdon  Non-Reflecting  Boundary  Conditions  to 
Shallow  Water  Models,  PhD  Thesis,  Naval  Postgraduate  School,  Monterey,  CA, 
2003. 


INITIAL  DISTRIBUTION  LIST 


1.  Defense  Technical  Information  Center 
Ft.  Belvoir,  Virginia 

2.  Dudley  Knox  Library 
Naval  Postgraduate  School 
Monterey,  California 

3.  Lt  Col  John  R.  Dea,  USAF 

Air  Force  Institute  of  Technology 
Wright-Patterson  AFB,  OH 

4.  CDR  Vince  J.  van  Joolen,  USN 
United  States  Naval  Academy 
Anapolis,  MD 

5.  Major  J.  M.  Lindquist,  USA 
Naval  Postgraduate  School 
Monterey,  CA 

6.  B.  Neta 

Naval  Postgraduate  School 
Monterey,  CA 

7.  F.  X.  Giraldo 

Naval  Postgraduate  School 
Monterey,  CA 


